Factorizing complex symmetric matrices with positive definite real and imaginary parts

نویسنده

  • Nicholas J. Higham
چکیده

Complex symmetric matrices whose real and imaginary parts are positive definite are shown to have a growth factor bounded by 2 for LU factorization. This result adds to the classes of matrix for which it is known to be safe not to pivot in LU factorization. Block LDLT factorization with the pivoting strategy of Bunch and Kaufman is also considered, and it is shown that for such matrices only 1× 1 pivots are used and the same growth factor bound of 2 holds, but that interchanges that destroy band structure may be made. The latter results hold whether the pivoting strategy uses the usual absolute value or the modification employed in LINPACK and LAPACK.

برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

منابع مشابه

On a Theorem of Haupt and Wirtinger concerning the Periods of a Differential of the First Kind, and a Related Topological Theorem

On a closed Riemann surface Ra of genus g there exist g linearly independent differentials of the first kind, wu • • • , ws,.and their integrals around 2g canonical cycles or retrosections, Oi, • • • , a„, bi, ••• , ba, are usually put together to form a gX2g matrix (an; ßn) = (A; B), i, j =1, • • • , g, where a,-,is the integral of w{ over a,. Riemann showed that the entries in this matrix wer...

متن کامل

No-gaps Delocalization for General Random Matrices

We prove that with high probability, every eigenvector of a random matrix is delocalized in the sense that any subset of its coordinates carries a non-negligible portion of its `2 norm. Our results pertain to a wide class of random matrices, including matrices with independent entries, symmetric and skew-symmetric matrices, as well as some other naturally arising ensembles. The matrices can be ...

متن کامل

On a General Solution of Partially Described Inverse Quadratic Eigenvalue Problems and Its Applications

In this paper, we consider to solve a general form of real and symmetric n× n matrices M , C, K with M being positive definite for an inverse quadratic eigenvalue problem (IQEP): Q(λ)x ≡ (λ2M + λC +K)x = 0, so that Q(λ) has a partially prescribed subset of k eigenvalues and eigenvectors (k ≤ n). Via appropriate choice of free variables in the general form of IQEP, for k = n: we solve (i) an IQE...

متن کامل

A FETI Method for a Class of Indefinite or Complex Second- or Fourth-Order Problems

The FETI-DP domain decomposition method is extended to address the iterative solution of a class of indefinite problems of the form (K−σM)x = b, and a class of complex problems of the form (K − σM + iσD)x = b, where K, M, and D are three real symmetric positive semi-definite matrices arising from the finite element discretization of either second-order elastodynamic problems or fourth-order pla...

متن کامل

An Effective Preconditioner for a PML System For Electromagnetic Scattering Problem

In this work we are concerned with an efficient numerical solution of a perfectly matched layer (PML) system for a Maxwell scattering problem. The PML system is discretized by the edge finite element method, resulting in a symmetric but indefinite complex algebraic system. When the real and imaginary parts are considered independently, the complex algebraic system can be further transformed int...

متن کامل

ذخیره در منابع من


  با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید

برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

عنوان ژورنال:
  • Math. Comput.

دوره 67  شماره 

صفحات  -

تاریخ انتشار 1998